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The  stabilisation  region  of  turbulent  non-premixed  flames  of  natural  gas  mixtures  burning  in  a  hot  and 
diluted  coflow  is  studied  by  recording  the  flame  luminescence  with  an  intensified  high-speed  camera. 
The  flame  base  is  found  to  behave  fundamentally  differently  from  that  of  a  conventional  lifted  jet  flame 
in  a  cold  air  coflow.  Whereas  the  latter  flame  has  a  sharp  interface  that  moves  up  and  down,  ignition  ker¬ 
nels  are  continuously  being  formed  in  the  jet-in-hot-coflow  flames,  growing  in  size  while  being  con¬ 
verted  downstream.  To  study  the  lift-off  height  effectively  given  these  highly  variable  flame  structures, 
a  new  definition  of  lift-off  height  is  introduced.  An  important  parameter  determining  lift-off  height  is 
the  mean  ignition  frequency  density  in  the  flame  stabilisation  region.  An  increase  in  coflow  temperature 
and  the  addition  of  small  quantities  of  higher  alkanes  both  increase  ignition  frequencies,  and  decrease  the 
distance  between  the  jet  exit  and  the  location  where  the  first  ignition  kernels  appear.  Both  mechanisms 
lower  the  lift-off  height.  An  increase  in  jet  Reynolds  number  initially  leads  to  a  significant  decrease  of  the 
location  where  ignition  first  occurs.  Higher  jet  Reynolds  numbers  (above  5000)  do  not  strongly  alter  the 
location  of  first  ignition  but  hamper  the  growth  of  flame  pockets  and  reduce  ignition  frequencies  in 
flames  with  lower  coflow  temperatures,  leading  to  larger  lift-off  heights. 
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1.  Introduction 

Turbulent  non-premixed  jet  flames  burning  in  a  hot  oxidiser 
stream  containing  combustion  products  (jet-in-hot-coflow  flames) 
are  relevant  to  clean  combustion  techniques  such  as  flameless 
combustion,  high-temperature  air  combustion  (HiTAC),  excess  en¬ 
thalpy  combustion  and  mild  combustion  [1-3].  These  techniques 
are  strongly  related  in  the  sense  that  they  all  recover  exhaust  gas 
heat  (for  instance  using  regenerators  or  recuperators)  and  rely  on 
high  recirculation  ratios  to  ensure  proper  mixing  of  one  or  both 
reactant  streams  with  the  flue  gases  (carrying  an  enthalpy  deficit) 
before  any  reaction  takes  place.  The  benefits  of  these  combustion 
techniques  are  the  inherent  flame  stabilisation  due  to  the  high  oxi¬ 
diser  temperature  and  the  flat  temperature  profiles  void  of  peaks, 
leading  to  low  NOx  emissions.  Furthermore,  efficiencies  are  high 
due  to  the  exhaust  gas  heat  recovery. 

The  industrial  and  environmental  relevance,  combined  with  the 
scientific  challenge  of  modelling  these  flames  is  the  reason  for  the 
amount  of  attention  this  type  of  combustion  has  recently  received 
in  both  detailed  experimental  [4-8]  and  numerical  [9-14]  studies 
on  laboratory-scale  setups. 

An  important  parameter  in  turbulent  non-premixed  flames  is 
the  lift-off  height.  In  conventional  jet  flames,  the  lift-off  height, 
or  flame  stabilisation  point,  is  the  average  axial  height  of  the  sharp 
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flame  interface.  For  jet-in-hot-coflow  flames,  different  criteria  to 
determine  lift-off  height  are  found  in  literature.  Usually,  a  certain 
threshold  level  for  an  averaged  quantity  is  defined,  for  instance 
temperature  [15],  OH  concentration  [4,16]  or  luminescence  [17]. 
As  gradients  of  average  quantities  in  jet-in-hot-coflow  flames  are 
generally  weak,  specific  choices  of  their  threshold  values  can  have 
great  impact  on  the  resulting  value  for  the  lift-off  height.  Moreover, 
determining  the  lift-off  height  from  time-averaged  fields  does  not 
reflect  the  complex  dynamics  of  the  lifted  flame. 

Several  theories  have  been  developed  to  explain  the  physical 
mechanisms  governing  lift-off  for  turbulent  jet  flames.  An  exten¬ 
sive  review  by  Lyons  [18]  divides  these  theories  and  concepts  in 
five  different  groups.  Summarising,  several  theories  rely  on  an 
equilibrium  between  flame  speed  and  flow  velocity,  the  flame 
speed  being  that  of  a  laminar  [19]  or  a  turbulent  [20]  premixed 
flame,  or  that  of  a  more  complex  2-D  flame  structure  [21],  while 
other  theories  emphasise  the  competition  between  chemical  and 
flow  time  scales  [22],  or  the  role  of  large-scale  flow  structures  in 
transporting  flame  leading  edges  [23].  This  list  does  not  include 
any  theory  that  is  specific  to  flameless  combustion.  In  the  recent 
study  of  lift-off  heights  that  specifically  included  flameless  com¬ 
bustion  [15],  a  generic  model  relating  flow  and  chemistry  time 
scales  was  used.  However,  it  has  been  pointed  out  in  several  stud¬ 
ies  [6,8,24,25]  that  the  physics  of  jet-in-hot-coflow  flames  is  rather 
different  from  that  of  a  conventional  lifted  jet  flame.  For  example, 
Gordon  et  al.  [8]  showed  by  means  of  joint  temperature,  OH  and 
CH20  imaging  the  distinctly  different  stages  of  flame  formation 
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Nomenclature 

<v 

spatial-  and  time  averaging  operators,  respectively 

tr 

the  value  of  b2  at  a  distance  and  time  lapse  from  an  iso- 

fa, 

function  describing  flame  presence  at  an  axial  height 

lated  ignition  event 

and  at  an  azimuthal  angle  in  time  (-) 

Vo 

mean  value  of  the  axial  velocity  component  with  which 

^2 

function  describing  flame  presence  at  an  axial  height  in 

flame  pockets  are  transported  (m/s) 

time  (-) 

Vbdick 

axial  velocity  of  the  trailing  edge  of  a  projected  flame 

d 

inner  diameter  of  the  fuel  jet  tube  (m) 

pocket  (m/s) 

/ign;l 

ignition  frequency  per  flame  surface  area  (m-2  s-1) 

V  front 

axial  velocity  of  the  leading  edge  of  a  projected  flame 

/ign;2 

ignition  frequency  per  axial  length  (m_1  s_1) 

pocket  (m/s) 

/ign 

ignition  frequency  (1/s) 

V]et 

bulk  velocity  of  the  jet  in  the  fuel  tube  (m/s) 

hi 

lift-off  height  based  on  b\  (m) 

^max;co 

maximum  velocity  of  the  coflow  (m/s) 

^2 

lift-off  height  based  on  b2  (m) 

Zbjnin 

the  axial  coordinate  where  ignition  kernels  first  occur  (m) 

/c, ,  k2 

non-dimensional  prefactors 

Rejet 

the  Reynolds  number  of  the  jet,  based  on  the  bulk  veloc- 

f  1  2 

distance  between  zb.min  and  h2  (m) 

ity  (-) 

Ptr 

probability  that  an  ignition  event  generates  a  burning 

T  CO 

temperature  of  the  coflow  (I<) 

flame  at  an  axial  distance  and  a  time  lapse  (-) 

7"max;co 

maximum  temperature  of  the  coflow  (I<) 

P„, 

time  averaged  value  of  b\  (-) 

2 

Taylor  length  scale  (m) 

P»2 

time  averaged  value  of  b2  (-) 

0 

derivative  of  the  expected  residence  time  of  the  axial 

r,z,(j> 

coordinates  of  the  cylindrical  coordinate  system 

projection  of  a  flame  pocket  with  respect  to  z  (s/m) 

in  methane  jet  flames  in  a  vitiated  coflow.  The  identified  sequence 
started  with  the  formation  of  CH20-rich  regions,  followed  by  the 
formation  of  small  (sub-millimetre)  OH-rich  kernels  that  grow  in 
size,  rise  in  temperature  and  eventually  lead  to  flame  structures. 
The  presence  of  ignition  kernels  has  been  demonstrated  in  other 
experiments,  for  instance  in  [24].  In  the  recent  review  by  Mastora- 
kos  [26]  an  overview  is  given  of  both  numerical  and  experimental 
work  in  the  field  of  autoigniting  flames.  The  combined  role  of 
chemistry  and  turbulence  is  discussed,  with  particular  attention 
to  the  role  of  strain  rates  at  mixture  fraction  levels  most  prone  to 
autoignition,  the  so-called  most  reactive  mixture  fraction. 

In  view  of  these  results,  there  appears  to  be  a  gap  between 
experimental  findings  on  the  physical  processes  that  are  responsi¬ 
ble  for  flame  stabilisation  in  jet-in-hot-coflow  flames  and  the  con¬ 
ventional  theories  regarding  lift-off  height  in  non-premixed 
turbulent  flames.  Therefore,  in  this  paper  an  attempt  will  be  made 
to  lay  a  theoretical  foundation  relating  the  lift-off  height  in  jet-in- 
hot-coflow  flames  to  the  observed  autoignition  and  convection 
processes.  An  accurate  definition  of  lift-off  height  will  be  formu¬ 
lated  first,  and  this  definition  will  be  used  to  compare  the  lift-off 
height  of  flames  in  the  Delft  Jet-in-Hot-Coflow  (DJHC)  burner  for 
a  range  of  parameters,  including  coflow  temperature,  jet  fuel  com¬ 
position  and  jet  Reynolds  number.  Secondly,  the  effect  of  these 
parameters  on  the  statistics  of  the  autoignition  process  will  be 
studied  in  detail.  The  relation  between  the  statistics  and  the  lift¬ 
off  height  is  then  used  to  explain  the  found  trends  for  lift-off 
height. 


ation  and  convective  heat  transfer)  the  flame  burns  in  circumstances 
resembling  those  in  a  furnace  in  flameless  combustion  conditions, 
with  peak  temperatures  in  the  flame  lower  than  those  of  a  natural 
gas  jet  flame  in  normal  air.  The  design  is  similar  to  the  Adelaide 
JHC  burner  [  5  ].  The  important  difference  is  that  the  secondary  burner 
is  not  fully  but  partially  premixed,  which  allows  for  seeding  of  the 
various  gas  flows  with  micron  size  particles  that  act  as  flow  tracers 
in  particle  image  velocimetry  and/or  laser  Doppler  anemometry.  In 
this  paper  a  z,  r-coordinate  system  will  be  used  with  its  origin  at 
the  centre  of  the  fuel  pipe  exit.  The  z-coordinate  is  measured  verti¬ 
cally  upward. 

2.2.  Case  description 

Both  cold-coflow  flames  (so-called  conventional  lifted  flames) 
and  hot  coflow  flames  have  been  studied.  The  cold-coflow  air  has 
a  temperature  of  293  K,  and  a  maximum  mean  velocity  vmaxico  of 
0.5  m/s.  Table  1  lists  the  characteristics  of  the  two  different  hot  co¬ 
flows  that  were  considered  in  this  study.  Tmax;co  represents  the 
maximum  temperature  in  the  coflow,  see  Section  3.1  for  further 
details.  The  coflow  temperature  in  the  DJHC-I  flame  is  about  80  K 
higher  than  that  of  the  DJHC-V  flame.  The  average  02  mass  fraction 
in  the  coflow  is  8.4%  in  the  case  of  DJHC-I  and  9.5%  in  the  case  of 
DJHC-V.  Three  different  fuels  were  used  in  the  central  jet:  Dutch 

A 

z 


2.  Experimental  setup 

2.2.  DJHC  burner 

Fig.  1  shows  a  sketch  of  the  upper  part  of  the  DJHC  burner.  This 
burner  is  designed  to  deliver  a  flame  that  mimics  the  important  char¬ 
acteristics  of  flameless  combustion.  It  creates  a  turbulent  diffusion 
flame  of  a  gaseous  fuel  in  a  coflowing  oxidiser  stream  of  high  temper¬ 
ature  carrying  little  oxygen.  A  fuel  jet,  emerging  from  a  long  (approx¬ 
imately  200  diameters)  tube  with  an  inner  diameter  d  of  4.5  mm, 
develops  in  a  coflow  of  hot  and  diluted  air.  This  coflow  is  generated 
by  a  secondary  burner  inside  an  annulus  with  a  diameter  of  82.8  mm. 
The  coflow  mixture  fraction  can  be  varied,  influencing  both  the  co¬ 
flow  temperature  and  oxygen  level.  Due  to  the  enthalpy  deficit  of 
the  coflow  (the  coflow  loses  heat  to  the  surroundings  through  radi- 
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Fig.  1.  Relevant  dimensions  of  the  upper  part  of  the  DJHC  burner.  Dimensions  are  in 
millimetres. 


E.  Oldenhof  et  al. / Combustion  and  Flame  157  (2010)  1167-1178 


1169 


Table  1 

Secondary  burner  flows  and  resulting  maximum  coflow  velocity  (measured  with  LDA) 
and  maximum  coflow  temperature  (measured  with  CARS)  of  the  studied  DJHC  flames. 


Case 

Fuel  (nl/min) 

Air  (nl/min) 

Tnax;co  (K) 

ZAnax;co  (m/ S) 

^02;co  (-) 

DJHC-I 

16.1 

224 

1540 

4.6 

8.4% 

DJHC-V 

15.3 

231 

1460 

4.3 

9.5% 

Table  2 

Compositions  of  the  three  fuels  used 

in  this  study. 

[Mole  %] 

Fuel  I 

Fuel  II 

Dutch  natural  gas 

n2 

15.0 

15.0 

14.4 

ch4 

85.0 

81.0 

81.3 

C2H6 

- 

4.0 

3.7 

Rest 

- 

- 

0.6 

natural  gas  and  two  synthetic  gases  resembling  Dutch  natural  gas, 
either  a  mixture  of  methane  and  nitrogen  (fuel  I)  or  a  mixture  of 
methane,  ethane  and  nitrogen  (fuel  II),  see  Table  2.  The  jet  Rey¬ 
nolds  number  was  varied  between  3000  and  9500,  and  is  based 
on  the  bulk  velocity  in  the  fuel  tube  and  the  dynamic  viscosity  of 
natural  gas  at  300  I<  for  the  cold-coflow  flames  (1.14  x  10  5  Pa  s) 
and  that  at  450  I<  for  the  hot-coflow  flames  (1.60  x  10  5  Pa  s).  Un¬ 
like  the  cold-coflow  flames,  the  hot-coflow  flames  do  not  have  to 
be  lit  externally. 

2.3.  Optical  setup  and  image  processing 

The  flame  structure  between z  =  0  mm  (corresponding  to  the  fuel 
jet  exit)  and  z  =  1 60  mm  was  studied  with  an  intensified  high-speed 
camera  (Lambert  Instruments  HI-CAM  CR).  This  camera  has  a  fibre- 
optically  coupled  two-stage  intensifier  with  a  sensitivity  exceeding 
50  mA/W  for  wavelengths  between  200  nm  and  550  nm.  It  was 
operated  at  a  gate  time  of  400  ps.  The  CMOS  sensor  has 
1280  x  1024  pixels,  but  only  a  quarter  (1280  x  256)  was  used  to 
achieve  higher  frame  rates.  A  Nikkor  UV  1 :4.5,/=  105  mm  lens  was 
used,  with  a  spectral  transmission  level  exceeding  70%  between 
220  nm  and  900  nm.  The  resolution  of  the  resulting  image  is  approx¬ 
imately?  pixels  mm-1.  The  field  of  view  was  wide  enough  to  capture 
the  width  of  the  flame  up  to  z  «  130  mm.  For  the  jet-in-hot-coflow 
flames,  1 04  frames  were  acquired  at  a  frame  rate  of  2000  fps  for  each 
case.  An  identical  frame  rate  and  gate  time,  albeit  at  a  lower  intensi¬ 
fier  gain,  was  used  to  study  the  dynamics  of  the  cold-coflow  flames. 
To  determine  the  lift-off  heights  of  the  cold-coflow  flames,  5000 
frames  were  taken  at  200  fps,  to  have  a  larger  averaging  time 
accounting  for  the  longer  time  scales  involved.  A  background  image 
was  subtracted  from  each  image.  This  background  image  was  con¬ 
structed  by  averaging  over  1000  frames  with  identical  exposure 
time,  frame  rate  and  intensification,  but  without  the  presence  of  a 
flame. 

To  process  the  images,  an  averaging  “disc”  with  a  radius  of  8  pix¬ 
els  was  applied  after  subtraction  of  the  background  image  to  reduce 
the  influence  of  noise.  Then  the  flame  boundaries  were  determined 
based  on  a  threshold  intensity  level  of 8  (255  being  the  maximum  va¬ 
lue),  making  use  of  the  built-in  Matlab  routine  bwboundaries.m.1 
The  dependence  of  the  resulting  lift-off  height  on  chosen  threshold 
values  is  small:  within  the  margin  of  reasonable  values  for  the  thresh¬ 
old  (flame  pockets  are  detected  as  such,  but  noise  is  rejected)  the  dif¬ 
ference  in  resulting  lift-off  is  3  mm. 

The  instantaneous  leading-  and  trailing  edge  speeds  of  flame 
pockets  vhack  and  ^front  are  determined  from  the  boundary  data  by 
matching  the  locations  of  the  projected  edges  in  sequential  images 


1  Supplementary  material  available,  showing  five  snapshots  before  image  smooth¬ 
ing  and  after  smoothing,  along  with  the  detected  boundaries. 


taken  at  tn  and  tn+ 1.  The  algorithm  that  determines  these  matches 
first  finds  candidates  for  the  front  edges  at  tn+ 1  from  those  at  tn, 
and  selects  the  one  with  the  smallest  positive  increment.  If  a  front 
edge  at  tn  has  multiple  matches  to  tn+ 1,  the  one  with  the  smallest 
increment  is  selected.  A  front  edge  with  no  forward  match  has 
merged  with  the  next  (downstream)  trailing  edge,  and  this  trailing 
edge  is  therefore  also  excluded  as  a  candidate  for  a  match  with  trail¬ 
ing  edges  at  tn+ 1  •  Front  edges  at  tn+ 1  that  have  no  backward  match  are 
recognised  as  being  formed  by  an  ignition  event,  and  both  the  front 
edge  and  the  trailing  edge  of  this  kernel  are  excluded  as  candidates 
for  a  match  with  edges  at  tn.  The  remaining  trailing  edges  at  tn  and 
tn+ 1  are  now  matched  in  order  of  their  axial  location.  As  a  result,  trail¬ 
ing  edges  can  have  both  positive  and  negative  speeds,  whereas  front 
edge  speeds  are  assumed  to  have  positive  speeds  always.  The  accu¬ 
racy  of  extracting  speeds  from  the  increments  is  limited  by  the  pre¬ 
cision  with  which  the  location  of  the  flame  boundary  can  be 
determined.  This  aspect  has  been  investigated  experimentally,  by 
projecting  a  light  dot  moving  on  a  screen,  recording  it  with  identical 
gain,  gate  time  and  frequency  and  processing  the  resulting  images 
with  the  boundary  detection  routine.  The  resulting  mean  speed 
was  found  to  be  unbiased,  as  it  approximates  the  true  speed  within 
1%.  The  95%  error  interval  of  individual  measurements  is  estimated 
at  ±1.1  m/s,  independent  of  the  light  dot  speed.  It  should  be  realised 
that  growth  speeds  of  ignition  kernels  at  the  moment  of  formation 
are  inherently  ill  defined,  because  of  the  sudden  increase  of  flame 
luminescence  over  their  extents. 

2.4.  LDA  system 

LDA  measurements  were  performed  with  a  two-component, 
dual  beam  TSI-system.  The  green  line  (514.5  nm)  and  blue  line 
(488  nm)  of  a  10  W  Continuum  Argon-ion  laser  were  used  to  mea¬ 
sure  the  axial  and  radial  velocity  components  directly.  Two  of  the 
incident  beams  (one  of  each  colour)  were  frequency  pre-shifted 
over  40  MHz  by  a  Bragg  cell  to  enable  the  detection  of  instanta¬ 
neous  flow  reversals.  The  length  and  diameter  of  the  measurement 
volume  were  1.7  mm  and  0.12  mm,  respectively.  Scattered  light 
was  collected  in  back-scatter  mode.  The  photomultiplier  output 
signals  were  electronically  down-mixed,  and  subsequently  col¬ 
lected  by  a  FSA-3000  signal  processor.  All  statistics  were  computed 
as  transit-time  weighted  results  to  eliminate  the  effects  of  the 
velocity  bias.  Autocorrelation  functions  of  the  axial  velocity  com¬ 
ponent  were  constructed  from  time  series  with  4  x  105  velocity 
samples  that  were  acquired  at  a  mean  data  rate  of  approximately 
500  Hz  by  using  the  slotting  method  with  local  normalisation  [27]. 

2.5.  CARS  system 

Temperatures  were  determined  with  a  CARS  system  that  has 
been  described  in  detail  in  [28].  In  a  planar-boxcars  phase-match¬ 
ing  configuration,  a  probe  volume  of  700  pm  length  and  35  pm 
diameter  is  obtained.  The  single-shot  imprecision  of  the  system 
(one  standard  deviation)  is  1-4%  over  a  range  from  2000  K  to 
300  K.  The  inaccuracy  (the  systematic  error)  is  estimated  to  be 
20  K.  For  each  position  in  space,  mean  temperatures  were  deter¬ 
mined  from  the  results  of  1000  single-shot  CARS  spectra. 

3.  Results  and  discussion 

3.1.  Details  of  the  coflow 

Fig.  2  shows  the  mean  temperature  field  at  z  =  3  mm  of  the  two 
coflows  DJHC-I  and  DJHC-V.  The  coflow  temperature  is  not  con¬ 
stant,  but  depends  on  the  radial  location  r. 
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Fig.  2.  Mean  temperature  in  the  DJHC-I  and  DJHC-V  flames  at  z  =  3  mm,  obtained 
with  coherent  anti-Stokes  Raman  spectroscopy  (CARS).  The  RMS  of  the  coflow 
temperature  is  typically  120  K.  Dutch  natural  gas  was  used  in  the  jet  stream  at  a  jet 
Reynolds  number  of  4500. 


Fig.  3.  Mean  axial  velocity  in  the  DJHC-I  and  DJHC-V  flames  at  z  =  3  mm,  obtained 
with  laser  Doppler  anemometry.  Dutch  natural  gas  was  used  in  the  jet  stream  at  a 
jet  Reynolds  number  of  4500. 


The  velocity  fields  of  the  DJHC-I  and  DJHC-V  flames  at  z  =  3  mm 
are  presented  in  Fig.  3.  The  larger  mass  flow  of  the  DJHC-V  coflow 
is  compensated  by  its  lower  temperature  and,  consequently,  higher 
density  to  yield  nearly  identical  velocities. 


3.2.  Visual  observations 

Fig.  4  shows  two  images  of  the  flame  base  of  a  conventional 
lifted  jet  flame  (left)  and  the  DJHC-V  (right).  In  both  cases  the  fuel 
is  Dutch  natural  gas  and  Rejet  =  4500.  Qualitative  differences  be¬ 
tween  the  conventional  lifted  flame  and  the  jet-in-hot-coflow 
flame  are  evident.  The  conventional  lifted  flame  has  a  sharp,  con¬ 
nected  interface  with  mild  indentations.  The  jet-in-hot-coflow 
flame  shows  large  spatial  and  temporal  variations  in  flame  inter¬ 
face  height.  Isolated  flame  pockets  develop  for  z  >  80  mm.  These 
grow  while  being  convected  downstream  where  they  merge  thus 
forming  a  more-or-less  continuous  flame  zone  at  z  «  120  mm. 

To  study  the  time  evolution  of  the  pockets,  the  boundary  detec¬ 
tion  routine  mentioned  earlier  will  be  used.  As  an  example,  Fig.  5 
shows  three  sequential  luminescence  images  (DJHC-V,  fuel  II, 
Rejet  =  7000)  with  the  green  contours  indicating  the  boundaries 


r  |  in  in  |  r  [min] 


Fig.  4.  Images  of  the  flame  base  of  a  conventional  lifted  flame  of  Dutch  natural  gas 
at  Rejet  =  4500  (left)  and  of  the  jet-in-hot-coflow  flame  DJHC-V  (Dutch  natural  gas 
at  Rejet  =  4500,  right). 


of  detected  flame  pockets.2  In  the  following,  regions  in  space  dis¬ 
playing  enough  luminescence  to  be  recognised  by  the  boundary 
detection  routine  as  burning  will  be  referred  to  as  flame  pockets. 
Flame  pockets  at  their  first  instant  of  detection  will  be  referred  to 
as  ignition  kernels.  The  grey  bars  on  the  right  of  each  image  in 
Fig.  5  represent  the  projection  of  the  flame  pockets  on  the  vertical 
coordinate  axis.  This  one-dimensional  function,  referred  to  as 
b2(z,t)  has  logical  values  (zero  if  no  flame  pocket  is  present  and 
one  if  a  flame  pocket  is  present).  Most  of  this  paper  will  be  related 
to  the  statistics  of  this  function  in  time.  Although  the  projection  on 
the  axis  leads  to  a  loss  of  information,  there  are  good  reasons  to 
study  b2  and  not  the  two-dimensional  boundary  information.  Most 
importantly,  due  to  the  axi-symmetry  of  the  flame,  the  relevant  igni¬ 
tion  statistics  are  a  function  of  z  only.  As  a  consequence,  b2(z,  t)  con¬ 
tains  all  information  needed  to  get  the  desired  information  on 
ignition  statistics.  Line-of-sight  effects  makes  the  distinction  be¬ 
tween  apparently  overlapping  pockets  impossible,  as  both  the  front 
and  the  back  side  of  the  flame  are  observed.  However,  statistics  on 
speeds  of  trailing  and  leading  edges  can  be  acquired  from  b2(z,t)t 
when  edges  of  flame  pockets  occur  axially  separated  from  other 
flame  pockets.  Therefore,  the  essential  features  of  the  observed  pro¬ 
cess,  namely  the  formation  of  ignition  kernels  and  the  growth  of  the 
resulting  flame  pockets,  can  be  retrieved  by  studying  b2(z,t)  alone. 

The  two  plots  in  Fig.  6  show  a  small  fraction  of  b2(z,  t)  for  a  jet- 
in-hot-coflow  flame  and  a  conventional  lifted  flame  (cold-coflow 
flame).  The  hot-coflow  flame  is  seen  to  behave  very  differently 
from  the  ordinary  jet  diffusion  flame,  in  which  a  single  sharp  inter¬ 
face  moves  up  and  down.  The  flame  pocket  behaviour  in  the  jet-in- 
hot-coflow  flames  bears  the  characteristics  of  a  process  of  random 
ignition  kernel  formation,  growth  and  convection. 

3.3.  Definition  of  lift-off  height 

The  lift-off  height  is  related  to  the  probability  of  the  presence  of 
flame  pockets.  Two  different  flame  probabilities  are  defined,  Pbl 
and  Pb2.  Pbi  (z)  is  the  probability  of  finding  a  flame  pocket  anywhere 


2  Supplementary  content  available,  showing  an  animation  of  the  flame  DJHC-V 
with  fuel  II  and  Rejet  =  7000  during  120  ms,  slowed  down  by  a  factor  of  400.  The  three 
frames  shown  in  Fig.  5  can  be  seen  at  t  =  6.0,  7.0  and  8.0  ms. 
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Fig.  5.  Three  sequential  images  from  the  case  DJHC-V,  fuel  II,  Rejet  =  7000,  showing  the  evolution  in  time  of  the  detected  flame  boundaries  (with  a  1  ms  time  separation, 
intermediate  frames  are  removed).  The  bar  on  the  right  shows  whether  a  flame  pocket  is  present  on  the  axial  location  (grey)  or  not  (white),  and  the  numbers  indicate  which 
trailing  edges  (numbers  preceded  by  a  “t”)  and  front  edges  are  matched  by  the  algorithm  that  determines  the  increments,  from  which  the  projected  flame  speeds  vhack  and 
front  are  calculated.  The  new  number  “82”  in  the  middle  picture  shows  that  an  ignition  event  was  recognised. 
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Fig.  6.  The  time  evolution  of  the  axial  location  of  flame  pockets  b2(z,  t).  (a):  DJHC-V, 
fuel  II,  Rejet  =  5000;  (b)  a  lifted  flame  in  cold  coflow,  Dutch  natural  gas,  Rejet  =  4500. 

on  a  radial  line  stretching  outward  from  the  burner  axis,  as  a  func¬ 
tion  of  axial  height.  It  could  be  determined  from  a  planar  imaging 
technique,  such  as  planar  laser-induced  fluorescence  (PLIF)  by 
observing  a  radial  cross-section  that  stretches  outward  from  the 
centre  and  calculating  the  fraction  of  time  a  burning  interface  is 
found  as  a  function  of  z.  PM  cannot  be  determined  directly  in  this 
experimental  setup,  since  the  luminescence  images  result  from  a 
line-integration.  However,  by  assuming  axi-symmetry  of  the  flame 
and  statistical  independence  of  the  “front”  and  “back”  of  the  flame, 
it  can  be  reconstructed  from  the  flame  boundary  observations  on 
the  centreline  of  the  image: 

PM(z)  =  l-(l-PM;d(z))1/2,  (1) 

where  Pb\;Ci{z)  is  the  fraction  of  time  a  pixel  at  height  z  on  the  centr¬ 
eline  contains  a  burning  pocket.  The  second  probability  Pb2(z)  is  that 


of  finding  a  flame  pocket  at  a  certain  axial  height,  and  is  the  expec¬ 
tation  of  the  function  b2(z,  t)  that  was  discussed  in  the  previous  sec¬ 
tion.  This  quantity  is  best  visualised  by  drawing  a  horizontal  line  in 
one  of  the  graphs  in  Fig.  6.  Pb2  is  the  fraction  of  the  line  that  is  in  the 
grey  region,  as  a  function  of  z. 

The  function  b2  is  equal  to  one  whenever  bi  (0,  z,  t)  is  equal  to 
one  for  an  azimuthal  angle  0,  and  it  is  as  such  not  an  azimuthally 
averaged  property.  However,  studying  its  time  evolution  has  a 
clear  advantage  over  studying  the  time  evolution  of  b\  since  it  is 
not  affected  by  “in-plane”  transport.  This  enables  a  more  rigorous 
analysis  of  the  ignition  and  transport  processes  and  for  this  reason, 
the  focus  will  be  primarily  on  b2.  In  this  respect,  the  used  technique 
has  an  advantage  over  a  planar  diagnostic  technique,  where  one 
cannot  with  certainty  distinguish  between  true  flame  islands  and 
filament-like  structures  being  convected  through  the  measure¬ 
ment  plane  [29],  and  simultaneous  measurement  of  the  velocity 
component  into  the  plane  is  needed  to  interpret  planar  images  of 
radicals  [30]. 

Both  probabilities  PM  and  Pb2  are  shown  in  Fig.  7  for  the  two 
flames  DJHC-I  and  DJHC-V.  The  lift-off  height  hi  is  now  defined 
as  the  axial  location,  where  Pbl  equals  50%.  A  second  lift-off  height 
h2  based  on  Pb2  will  also  be  used  in  the  analysis.  The  locations  cor¬ 
responding  to  these  lift-off  heights  are  shown  together  with  aver¬ 
age  luminescence  images  in  Fig.  8.  As  expected,  hi  is  located  at  the 
height  where  the  gradient  in  luminescence  is  high.  Note  that  in  the 
two  flames  at  the  bottom  of  the  picture,  there  is  a  large  margin  for 
choosing  a  lift-off  height  based  on  time  averaged  visual 
observations. 

Figs.  9  and  10  show  the  lift-off  heights  as  a  function  of  the  jet 
Reynolds  number  for  the  different  flames.  Both  definitions  for 
lift-off  height  yield  similar  trends.  The  minimum  lift-off  height 
hi  =  18d  is  found  in  flame  DJHC-I  with  natural  gas  as  fuel  and  at 
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Fig.  7.  The  probabilities  of  finding  a  flame  at  an  axial  height  on  any  location  (Pb2) 
and  at  an  axial  height  on  a  radial  line  (Pm),  for  the  DJHC-I  and  DJHC-V  flames  with 
natural  gas,  Rejet  =  4500.  The  symbols  indicate  the  location,  where  the  probabilities 
are  0.5. 


Fig.  9.  Lift-off  heights  based  on  Pfal  against  jet  Reynolds  number.  Note  that  the  lift¬ 
off  height  of  case  DJHC-V  with  fuel  I  exceeded  160  mm  for  each  jet  Reynolds 
number.  The  thick,  dashed  line  represents  the  correlation  for  lift-off  height  by 
Kalghatgi  [20]. 
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Fig.  8.  The  locations  of  Pbl  =  50%  (hi)  and  of  Pb2  =  50%  ( h2 )  marked  on  averaged 
images  of  four  flames. 


a  Reynolds  number  of  4500.  In  all  cases,  the  value  of  h\  (or  h2 )  has  a 
minimum  at  a  moderate  Reynolds  number  of  approximately  5000. 


The  value  of  h2  exceeds  the  field  of  view  of  the  camera  for  the  flame 
DJHC-V  with  fuel  I  at  Rejet  =  9500.  This  happens  in  several  more 
cases  for  hi,  most  notably  in  flame  DJHC-V  with  fuel  I  at  all  Rey¬ 
nolds  numbers.  The  difference  between  h2  and  hi  for  identical 
cases  increases  with  larger  Reynolds  numbers.  This  trend  will  be 
explained  in  Section  3.8. 

To  place  these  results  in  some  context,  two  reference  cases 
should  be  mentioned.  The  lift-off  height  in  the  vitiated  coflow  bur¬ 
ner  as  described  in  Cabra  et  al.  [4]  has  been  studied  for  a  range  of 
Reynolds  numbers,  with  a  CH4/air  fuel  stream  and  Tco  equal  to 
1350  K  [17],  and  with  an  H2/N2  fuel  stream  and  Tco  equal  to 
1045  K  [12].  The  fuel  jet  has  an  inner  diameter  d  of  4.57  mm,  com¬ 
parable  to  the  4.5  mm  of  the  DJHC  burner.  In  the  first  study,  lift-off 
heights  (based  on  average  flame  luminescence)  at  a  coflow  velocity 
of  4.2  m/s  ranged  between  h/d  «  11  to  h/d  ^  30  for  a  jet  Reynolds 
number  ranging  between  23,000  and  72,000  ( v]et  from  80  m/s  to 
250  m/s),  in  a  linear  fashion.  In  the  second  study,  lift-off  heights 
ranged  from  h/d  ^  3  to  h/d  «  20  for  a  Reynolds  number  range  of 
1 1 ,000  to  35,000  ( Pjet  from  45  m/s  to  1 60  m/s)  with  a  coflow  veloc¬ 
ity  equal  to  3.5  m/s.  A  monotonically  increasing  trend  was  found 
again,  with  a  steepening  towards  higher  jet  velocities.  An  initially 
decreasing  trend  as  seen  here  has  not  been  found  in  either  case, 
the  coflow  temperature  profile  in  these  cases  is  however  flat  and 
the  Reynolds  numbers  involved  are  higher. 

Variation  in  fuel  composition  through  the  addition  of  higher  al¬ 
kanes  (ethane  in  fuel  II  and  propane  and  butane  in  Dutch  natural 
gas)  affects  lift-off  height  strongly.  This  is  in  agreement  with  the 
strong  influence  of  chemistry  in  jet-in-hot-coflow  flames,  which 
has  been  pointed  out  in  previous  numerical  studies  [10,31  ].  The  re¬ 
sults  are  in  line  with  the  fact  that  higher  alkanes  tend  to  reduce  the 
autoignition  delay  times  of  natural  gas  mixtures  [32]. 

Figs.  9  and  10  also  show  the  lift-off  behaviour  of  conventional 
jet  flames  for  the  three  different  fuels.  In  this  case,  the  dependence 
of  lift-off  height  on  fuel  composition  is  very  weak.  The  Reynolds 
number  dependence  as  indicated  with  the  thick  dashed  line  fol¬ 
lows  the  relation  of  Kalghatgi  [20]  qualitatively,  being  linearly 
dependent  on  the  Reynolds  number.  The  coflow  velocity  of 
0.5  m/s  is  accounted  for  in  this  calculation  through  the  effective 
velocity  as  proposed  by  Montgomery  et  al.  [33]. 

The  issue  to  be  addressed  now  is  what  processes  cause  the 
strong  differences  in  lift-off  height  for  different  fuels,  Reynolds 
numbers  and  coflow  temperatures.  For  this  purpose,  the  probabil¬ 
ity  function  Pb2(z)  will  be  studied. 
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Fig.  10.  Lift-off  heights  based  on  Pb2  against  jet  Reynolds  number.  See  for  the  Fig.  12.  The  location  of  first  ignition  events  zb.min  against  jet  Reynolds  number, 

legend  Fig.  9. 


3.4.  Axial  location  of  first  occurrence  of  autoignition  kernels 

Fig.  11  illustrates  the  influence  of  the  jet  Reynolds  number  on 
the  probability  Pb2  and  thus  on  the  lift-off  height  h2.  It  can  be  seen 
that  the  decrease  of  lift-off  height  for  increasing  Reynolds  number 
up  to  Rejet  =  5000  is  mainly  due  to  a  horizontal  shift  of  the  curve. 

The  location  where  ignition  kernels  first  occur,  zb,min  (arbitrarily 
defined  as  the  location,  where  Pb2  equals  2.5  x  1(T3),  shifts  down¬ 
ward  and  so  does  the  entire  curve.  For  higher  Reynolds  numbers, 
Zb;min  IS  hardly  affected,  but  the  gradient  of  the  curve  of  Pb2 
decreases. 

Fig.  12  illustrates  the  effect  of  coflow  temperature,  fuel  compo¬ 
sition  and  Reynolds  number  on  zb,min.  Replacing  only  4%  methane 
by  ethane  (fuel  II  vs.  fuel  I)  strongly  decreases  zfa;mi-n,  especially  in 
DJHC-V  where  the  coflow  temperature  is  relatively  low.  The  impact 
of  increasing  the  jet  Reynolds  number  is  remarkable.  Instead  of 
transporting  gas  mixtures  further  away  from  the  jet  exit  before 
autoignition  due  to  the  higher  velocities,  it  causes  zb.min  to  decrease 
initially  in  all  cases.  In  the  flame  DJHC-V,  with  fuel  I,  zb.min  rises 
again  above  Rejet  =  5000,  indicative  of  hindered  formation  of  igni¬ 
tion  kernels.  Simply  translating  the  axial  coordinate  to  the  resi¬ 
dence  time  of  a  fluid  parcel  leads  to  the  conclusion  that 
autoignition  delay  times  are  reduced  under  the  influence  of 
increasing  jet  Reynolds  number.  Under  the  condition  that  the  igni- 


Fig.  11.  The  curve  of  Pb2  for  a  range  of  Reynolds  numbers,  DJHC-V  fuel  II. 


tion  time  is  larger  than  the  typical  turbulence  time,  faster  mixing 
will  lead  to  earlier  generation  of  well-mixed  spots  promoting 
autoignition,  a  possibility  discussed  in  [26].  This  conclusion  can 
however  not  be  drawn  here,  because  turbulence  will  cause  varia¬ 
tion  in  fluid  parcel  residence  times  at  a  given  location.  Moreover, 
the  faster  entrainment  of  the  hotter  inner  part  of  the  coflow  for 
increasing  Reynolds  number  is  expected  to  play  a  dominant  role. 

3.5.  Transport  of  flame  pockets 

After  a  flame  pocket  has  been  formed,  the  pocket  will  grow 
while  being  convected  downstream.  Clearly,  the  axial  velocity  of 
the  leading  edge  of  the  flame  pocket  vfront  must  be  larger  than 
the  axial  velocity  of  the  trailing  edge,  vback,  for  the  axial  extent  of 
the  flame  pocket  to  increase.  If  flame  stabilisation  by  flame  propa¬ 
gation  does  not  play  a  role,  vback  should  be  larger  than  zero.  The 
physical  relevance  of  these  velocities  to  the  lift-off  height  can  be 
demonstrated  in  a  simple  example  where  ignition  events  take 
place  at  a  fixed  location  and  at  regular  time  intervals  with  fre¬ 
quency  /ign.  At  streamwise  distances  from  the  ignition  location  lar¬ 
ger  than 

A z=(— - — )  (2) 

\  ^back  Vfront  /  /ign 

there  will  be  a  flame  at  any  instant  in  time.  This  illustrates  that 
increasing  vback  increases  the  lift-off  height,  while  increasing  the 
ignition  frequency  fign  decreases  the  lift-off  height. 

Histograms  showing  the  measured  values  of  vhack  and  ^front  of  a 
DJHC-V  flame  and  the  trailing  edge  speeds  of  a  conventional  lifted 
flame  determined  by  the  routine  as  described  in  Section  2.3  are 
shown  in  Fig.  13.  These  flame  speeds  are  not  to  be  confused  with 
the  propagation  speeds  of  flame  structures  such  as  studied  in  flat 
flame  configurations:  firstly  the  speed  is  not  measured  relative  to 
its  surrounding  fluid,  and  secondly  it  is  subject  to  horizontal  pro¬ 
jection.  The  correlation  between  these  speeds  and  axial  height  z 
is  very  weak.  As  a  consequence,  the  mean  values  of  the  pocket  edge 
speeds  can  be  treated  to  a  first  approximation  as  being  indepen¬ 
dent  of  axial  height  z. 

Attention  will  now  be  given  to  the  influence  of  the  different 
parameters  on  the  mean  axial  flame  speeds  /’front  and  vback.  The  re¬ 
sults  for  all  cases  are  shown  in  Fig.  14.  The  effect  of  jet  Reynolds 
number  on  /’front  is  small,  whereas  vback  clearly  increases.  This 
can  be  interpreted  as  a  mild  increase  of  the  mean  speed  with  which 
flame  pockets  are  convected  v0,  (assumed  to  be  the  average  of  vhack 
and  /’front),  combined  with  a  stronger  decrease  in  the  mean  pro- 
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Fig.  13.  Normalised  2-D  histograms  of  velocity  and  location  (greyscale  plot)  and 
pdfs  of  velocity  (red  line,  right  axis),  of  trailing  edge  in  a  conventional  lifted  flame 
(Dutch  natural  gas)  (a),  leading  edges  (b)  and  trailing  edges  (c)  of  flame  pockets  in 
flame  DJHC-V,  fuel  II.  The  jet  Reynolds  number  is  5000  in  all  cases.  The  2-D 
histogram  is  normalised  such  that  the  cell  with  the  maximum  number  of  counts  for 
each  case  gets  a  value  equal  to  one. 


jected  flame  propagation  speed  (the  difference  between  Pback  and 
^front  divided  by  two),  for  increasing  jet  velocity.  The  first  trend  is 
in  agreement  with  the  mean  location  of  the  reaction  zone  being 
near  the  outer  edge  of  the  jet.  CARS  measurements  indicate  that 
the  reaction  zone  at  z=120  mm  is  located  between  r=13mm 


Fig.  14.  Observed  front-  and  trailing  edge  speeds  in  axial  direction  at  different  jet 
Reynolds  numbers.  Both  the  tVck-  and  z>front  errorbars  are  based  on  the  RMS  value  of 
differences  between  the  values  obtained  from  the  first  and  last  set  of  5000  images 
of  their  respective  groups.  Cases  with  insufficient  statistics  (less  than  500  samples) 
were  discarded. 


and  r=  15  mm  (case  DJHC-I,  Rejet  =  4500  and  Dutch  natural  gas 
as  fuel).  Here,  the  mean  velocity  is  mainly  set  by  the  coflow  and 
is  thus  only  mildly  dependent  on  the  jet  Reynolds  number.  For 
example,  the  measured  mean  velocity  at  z=120  mm  and 
r=14mm  increases  from  5.0  m/s  at  a  jet  Reynolds  number  of 
3000,  to  6.6  m/s  at  Rejet  =  9500,  close  to  the  values  of  v0.  The  loca¬ 
tion  of  the  reaction  zone  also  explains  the  lack  of  correlation  be¬ 
tween  axial  height  and  flame  speed,  as  the  axial  gradients  of 
mean  velocity  in  the  coflow  are  small.  The  properties  of  the  flow 
field  in  the  flame  region  are  shown  in  more  detail  in  Section  3.7. 
The  DJHC-V  flames  have  consistently  lower  projected  flame  propa¬ 
gation  speeds,  which  might  be  caused  by  lower  true  flame  propa¬ 
gation  speeds  or  by  an  increase  of  curvature  of  the  contours 
along  which  the  flame  propagates.  The  latter  might  be  caused  by 
the  higher  oxygen  content  in  the  coflow,  increasing  the  stoichiom¬ 
etric  mixture  fraction  and  thereby  moving  the  reaction  zone  closer 
to  the  shear  layer  of  the  jet. 


3.6.  Flame  transfer  probability 


As  explained  in  the  previous  section,  the  flame  pocket  growth  is 
determined  by  leading  and  trailing  edge  speeds.  The  function  that 
completely  describes  the  convection,  growth  and  possible  extinc¬ 
tion  of  newly  formed  ignition  kernels  is  the  flame  transfer  proba¬ 
bility,  Ptr(Az,At).  This  function  describes  the  probability  that  an 
ignition  event  J(z,  t )  generates  a  burning  flame  at  an  axial  distance 
A z  and  time  lapse  At  from  this  event.  In  the  case  of  fixed  flame 
speeds, 


^fr;fixed  (^Z,  At) 


1  if  ^back  ^  At  ^  ^front 

0  otherwise. 


This  expression  for  Ptr  ignores  the  possibility  of  velocity  variations, 
or  extinction,  of  a  newly  formed  flame  pocket.  Therefore,  a  more  re¬ 
fined  method  is  needed,  using  b2(z,t)  directly  to  obtain  the  flame 
transfer  probability. 

To  demonstrate  the  relation  between  b2  and  Ptr,  consider  a  sin¬ 
gle  realisation  of  b2(z  +  Az,  t  +  At),  conditional  on  an  ignition  event 
J(z,t).  The  resulting  field  of  b2(z  + Az,  t  + At)  is  the  union  of  the 
field  in  absence  of  this  ignition  with  the  flame  evolution  of  this 
specific  ignition  event,  tr(Az,  At).  The  expectation  of  b2  (the  argu¬ 
ments  are  omitted  for  brevity)  conditional  on  the  ignition  event  I 
follows  from 

E{b2|I}  =  E{tr }  +  (1  -  E{tr})E{62}  =  Ptr  +  (1  -  Pfr)Pb2,  (4) 

using  statistical  independence  of  the  individual  ignition  events. 
Note  that  E  is  the  “expected  value”  operator  and  “|”  denotes  that 
a  conditional  expectation  is  considered.  Rewriting  Eq.  (4)  leads  to 
the  following  equation  (written  out  fully): 

d  ( k-t  AM_E{b2(z  +  Az,t  +  At)|I(z,t)}-Pij2(z  +  Az)  ^ 

MA2’At)~  1-Pt2(z  +  Az)  '  (5) 

In  practice  this  means  extracting  samples  of  b2(z,  t)  (see  Fig.  6) 
around  a  large  number  of  ignition  events  accounting  each  time 
for  the  locally  expected  value  of  b2  which  is  Pb2(z  + Az).  This  ap¬ 
proach  works  well,  provided  that  the  value  of  Pb2  is  not  too  large, 
and  ignition  events  can  be  identified  accurately.  Both  requirements 
are  fulfilled  by  considering  only  ignition  events  at  locations  of  low 
flame  probability,  say  Pb2  <  0.03.  The  resulting  statistics  consist 
typically  of  around  200  samples,  with  a  minimum  of  123  samples 
for  the  case  DJHC-V  with  fuel  I  at  a  jet  Reynolds  number  of  7000. 
Fig.  15  shows  the  flame  transfer  probability  Ptr(Az,  At)  for  four  Rey¬ 
nolds  numbers  in  DJHC-V  flames  with  fuel  II.  These  figures  clearly 
show  the  mean  downstream  convection  and  growth  of  newly  cre¬ 
ated  flame  pockets.  The  dotted  blue  lines  according  to  Eq.  (3)  corre¬ 
spond  well  to  the  found  flame  transfer  probability  functions.  For  the 
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Fig.  15.  The  flame  transfer  probability  PCn  for  flame  DJHC-V,  fuel  II,  with  Reynolds 
number  3000  (a),  5000  (b),  7000  (c)  and  9500  (d).  It  denotes  the  ensemble-averaged 
evolution  of  a  flame  kernel  starting  from  its  creation  at  At  =  0.  Each  plot  consists  of 
approximately  200  of  these  samples.  Note  that  at  the  highest  Reynolds  number,  Ptr 
becomes  more  sparse,  which  is  caused  by  occasional  extinction  events.  The  dotted 
blue  lines  show  the  pocket  evolution  corresponding  to  the  values  of  z>front  (steepest 
curve)  and  T’back.  as  determined  in  Section  3.5. 


highest  Reynolds  numbers,  the  effects  of  inhibited  growth  of  igni¬ 
tion  kernels  and  even  extinction  become  important,  as  can  be  ob¬ 
served  by  the  decrease  of  the  transfer  probability  in  the  region 
between  the  blue  dotted  lines  for  increasing  values  of  A z  and  At. 

The  function  Ptr  acts  as  a  convolution  kernel  on  the  ignition  fre¬ 
quency  density /ign;2  (2)  (the  frequency  of  formation  of  ignition  ker¬ 
nels  per  unit  axial  length  in  absence  of  a  flame)  to  yield  the  flame 
probability  Pb2  (with  the  introduction  of  the  dummy  variables  z! 
and  t'): 


Pb 2  (z,  t)  =  b2  =  1  -  exp  [-/ign;2  (z',  0  *  Ptr(Az,  At)] ,  (6) 

where  is  the  convolution  operator.  Therefore,  the  results  in 
Fig.  1 5  enable  the  quantification  of  the  growth  of  kernels  and  to  re¬ 
late  ignition  frequencies  to  the  flame  probability  Pb2. 

Using  the  time  invariance  of  the  ignition  frequency,  and  the 
ansatz 


J  Ptr(Az,  At)  dt  =  0A z, 


4 


1 


0  2  4  6  8  10  12  14  16 

Az  [mm] 


- 1 - 1 - 1 — 

-A —  DJHC-I,  fuel  I,  Re=5000 
—  DJHC-I,  fuel  II,  Re=5000 
-■ —  DJHC-V,  fuel  II,  Re=3000 
-• —  DJHC-V,  fuel  II,  Re=5000 
^ —  DJHC-V,  fuel  II,  Re=7000 
-x —  DJHC-V,  fuel  II,  Re=9500 


— i - 1 - 1 - r 

—  —  —  d^fjxed  A  Z  +  C 

Re=3,000 
Re=9,500 


Fig.  16.  The  value  of  the  time  integral  over  Ptr  as  function  of  Az,  for  several  cases. 
The  thick  dashed  black  lines  denote  the  slope  based  on  Eqs.  (7)  and  (9)  and  the 
flame  speeds  z>fr0nt  and  TWk,  of  flame  DJHC-V  with  fuel  II  for  a  jet  Reynolds  number 
of  3000  (top  line)  and  9500  (bottom  line).  The  thin  grey  lines  represent  the  values  of 
all  hot  coflow  cases,  to  give  an  impression  of  the  amount  of  scatter  between  cases. 


nolds  numbers  are  comparable  to  those  calculated  assuming  fixed 
flame  speeds,  in  which  case  0  simply  follows  from: 


^fixed 


1  1 

Pback  Pfront 


A  characteristic  of  all  curves  is  that  they  do  not  intersect  the  origin. 
This  is  caused  by  the  fact  that  ignition  kernels  already  have  a  certain 
size  at  the  moment  they  are  detected  by  the  image  processing  rou¬ 
tine.  They  appear  as  vague  “streaks”  of  a  certain  size  one  or  two 
frames  before  fully  igniting  and  being  accurately  detectable.  The 
average  kernel  length  in  axial  direction  upon  detection  is  between 
3  and  5  mm.  This  is  comparable  to  the  radius  of  the  averaging  disc 
(1  mm)  and  the  convective  displacement  during  the  400  ps  expo¬ 
sure  time  (2.4  mm).  Quantification  of  the  true  size  of  kernels  at 
the  moment  of  detection  is  therefore  not  undertaken  here. 


3.7.  Ignition  frequencies 

The  ignition  frequency  may  be  determined  from  the  flame 
boundary  data.  However,  validation  of  the  used  ignition  kernel 
detection  routine  against  Monte  Carlo  simulations  showed  that  it 
was  unreliable  in  regions,  where  Pb2  has  a  significant  value,  say 
Pb2  >  0.3,  and  at  higher  ignition  frequency  densities.  Another 
shortcoming  of  the  detection  routine  is  that  it  is  unable  to  recog¬ 
nise  ignition  events  in  rapid  succession  over  small  distances.  An 
exact  expression  for  /jgn;2(z)  derived  from  Eq.  (8)  involves  the  sec¬ 
ond  spatial  derivative  of  Pb2l  which  results  in  very  noisy  data. 
The  following  measure  for  the  mean  ignition  frequency  is  therefore 
used  instead,  based  on  the  more  robust  quantity  dPb2/dz.  Differen¬ 
tiating  Eq.  (8)  gives 


Eq.  (6)  can  be  rewritten  as  a  simple  integral  over  the  axial  height  z: 


Pb;2(z)  =  1  -  exp 


•<*>  J  /ign;2(z')(Z-Z')dz' 


The  parameter  0  is  the  derivative  of  the  expected  residence  time  of 
a  projected  flame  pocket  pertaining  to  an  ignition  event  (the  time 
between  arrival  of  the  front  edge  and  arrival  of  the  trailing  edge) 
with  respect  to  z.  Note  that  0  =  0  implies  no  flame  growth  and 
hence  no  flame  stabilisation  whilst  0  — >  00  corresponds  to  stabilisa¬ 
tion  by  flame  propagation.  Fig.  16  shows  the  time  integral  of  Ptr  for  a 
range  of  distances  from  the  ignition  event,  Az.  The  assumption  of 
linearity  in  Az  appears  to  be  appropriate.  The  values  for  lower  Rey- 


Pb2=l/2 


0 

2 


<Pb2= 1/2) 


/ign;2  dz 


0 

2"  (^2  —  Zb;min)  (/ign;2 )  • 


Renaming  the  length  scale  (h2  -  zb:min)  to  Li/2  gives: 


2 

dPb2 

£  1/2 

dz 

Pb2=V2 


(11) 


Fig.  17  shows  the  mean  ignition  frequency  (fign;2)  as  a  function  of  jet 
Reynolds  number  for  several  flames.  The  value  of  (/ign;2)  has  been 
determined  with  Eq.  (11)  using  linear  fits  to  the  curves  of  Pb2(z) 
and  those  in  Fig.  16  to  obtain  their  derivatives.  The  ignition  frequen¬ 
cies  thus  found  range  between  2  x  103  nrr1  s_1  and  6  x  104  m_1s_1. 
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Fig.  17.  Mean  ignition  frequency  densities  (fjgn;2)  as  determined  with  Eq.  (11) 
against  jet  Reynolds  number,  presented  on  a  semi-log  scale. 


The  average  ignition  frequency  density  calculated  by  the  image  pro¬ 
cessing  routine  reproduces  the  trends  and  figures  for  the  DJHC-V 
flames,  but  the  higher  values  are  not  reproduced,  which  is  due  to 
the  shortcomings  mentioned  earlier.  The  mean  ignition  frequency 
densities  of  flames  DJHC-I  are  at  least  two  times  larger  than  those 
in  flames  DJHC-V,  and  fuels  with  higher  alkanes  yield  consistently 
higher  values.  The  jet  Reynolds  number  hardly  has  any  influence 
on  the  ignition  frequencies  in  flame  DJHC-I,  whereas  flame  DJHC- 
V  is  more  sensitive  to  the  Reynolds  number,  especially  for  fuel  I. 

The  diminished  frequency  of  ignition  events  (in  DJHC-V)  and  in¬ 
creased  probability  of  kernel  extinction  and  limited  flame  pocket 
growth  (in  all  cases)  is  expected  to  originate  from  the  larger  strains 
in  the  turbulent  field  at  higher  Reynolds  numbers.  The  large  role  of 
strains  on  autoignition  has  been  pointed  out  in  previous  DNS-stud- 
ies  [34-36].  Echekki  et  al.  [36]  showed  the  evolution  of  ignition 
kernels  in  2-D  DNS  calculations  with  detailed  hydrogen-air  chem¬ 
istry.  A  species-specific  Damkohler  number  was  introduced,  and 
observed  extinction  events  were  related  to  excessive  heat  and  spe¬ 
cies  dissipation.  Very  recently,  experimental  evidence  has  been  ob¬ 
tained  showing  that  ignition  kernels  with  super-equilibrium  OH- 
levels  strongly  favour  locations  with  low  temperature  gradients 
[37].  In  order  to  verify  that  the  lower  ignition  frequencies  indeed 
originate  from  higher  levels  of  scalar  dissipation,  detailed  statistics 
on  the  history  of  the  scalar  field  that  ignition  kernels  experience 
prior  to  ignition  would  be  needed,  an  experimental  effort  not 
undertaken  in  this  study.  However,  some  insight  on  the  influence 
of  the  jet  Reynolds  number  on  the  turbulent  field  in  the  reaction 
zone  would  be  useful.  Therefore,  the  quantity  ( du'z/dz )2  was  deter¬ 
mined  from  LDA  measurements  by  first  computing  the  autocorre¬ 
lation  function  (acf)  of  the  axial  velocity  fluctuations  and  then 
invoking  Taylor’s  hypothesis  to  convert  time  delay  into  spatial  sep¬ 
aration  [38].  The  resulting  acfs  that  were  measured  at  several  ra¬ 
dial  locations  in  the  DJHC-I  flame  at  z  =  120  mm  for  Rejet  =  4500 
and  Rejet  =  9500  are  shown  in  Fig.  18. 

The  intercept  of  the  parabola  with  the  Az-axis  represents  the 
Taylor  length  scale  2  which  is  related  to  the  strain  rate  as  in 


The  resulting  values  for  the  Taylor  length  scale  are  listed  in  Table  3. 
The  Taylor  length  scale  decreases  (and  the  variance  of  the  velocity 
fluctuations  increases)  with  increasing  Reynolds  number.  This  is 
clear  evidence  that  the  mean  strain  rate  in  the  vicinity  of  the  reac¬ 
tion  zone  is  much  larger  at  the  higher  jet  Reynolds  number.  Both 
the  impact  of  jet  Reynolds  numbers  on  ignition  frequency  densities 


Fig.  18.  Autocorrelation  functions,  DJHC-I,  Rejet  =  4500  and  Rejet  =  9500.  The 
connected  and  dotted  line  represent  the  parabolic  fit  through  the  autocorrelation 
function  at  Az  =  0  for  Rejet  =  4500  (r  =  14  mm)  and  Rejet  =  9500  (r  =  15  mm), 
respectively.  These  are  the  locations  where  the  flame  resides  on  average. 


Table  3 

Flow  characteristics  in  DJHC-I  flames  at  three  radial  locations  for  Rejet  =  4500  (Dutch 
natural  gas)  and  Rejet  =  9500  (fuel  II).  The  Taylor  length  scales  between  parentheses 
should  be  addressed  with  some  care  because  of  the  high  turbulent  intensity  (above 
30%). 


Rejet  =  4500  Rejet  =  9500 


r  (mm) 

13 

14 

15 

13 

14 

15 

u  (m/s) 

5.95 

5.42 

5.31 

7.65 

6.57 

5.99 

\/ ii'2  /Ti  f-'l 

18.6% 

14.8% 

12.7% 

31.3% 

30.6% 

26.4% 

1  (mm)  4.8  [±0.2]  5.4  [±0.3]  5.9  [±0.3]  (3.8  [±0.2])  (3.4  [±0.2])  3.4  [±0.2] 


in  the  DJHC-V  flames  as  shown  in  Fig.  1 7  and  the  hindered  kernel 
growth  as  was  evident  in  Fig.  1 5  can  therefore  be  linked  to  changes 
in  the  turbulent  flow  field. 


3.8.  Evaluation  and  implications  on  lift-off  height 

In  the  analysis  of  the  mechanisms  governing  lift-off,  the  func¬ 
tion  b2(z,t)  was  considered  to  be  the  result  of  a  simple  stochastic 
process,  namely  that  of  randomly  occurring  local  events  generating 
infinitesimal  ignition  kernels  with  an  expected  flame  growth  that 
is  linear  in  time,  determined  from  the  flame  transfer  probability 
function  Ptr.  This  probability  Ptr  was  assumed  to  be  independent 
of  the  axial  location  z  in  the  stabilisation  region.  Several  deviations 
from  the  idealised  process  can  be  witnessed  however.  Ignition  ker¬ 
nels  do  not  appear  as  points,  but  as  small  regions  in  space,  with  a 
fast,  but  not  instantaneous,  transition  from  a  non-burning  to  a 
burning  state.  Ignition  events  are  unlikely  to  be  statistically  inde¬ 
pendent,  as  occasionally  successive  ignition  events  are  witnessed 
in  a  certain  region. 

Although  these  deviations  will  certainly  have  an  influence,  the 
essential  features  of  the  flames  (for  instance  as  shown  in  Fig.  15) 
strongly  resemble  those  of  the  proposed  process.  Stabilisation  by 
flame  propagation,  which  forms  the  basis  of  most  theories  on 
lift-off  in  turbulent  non-premixed  flames,  as  the  dominant  mecha¬ 
nism  is  out  of  the  question,  because  the  flame  speed  at  the  trailing 
edges  is  insufficient  to  maintain  the  flame  base  at  a  fixed  point  (see 
Fig.  14).  This  is  not  to  say  that  the  flame  speed  is  irrelevant,  since  it 
affects  the  growth  of  flame  pockets  and  thereby  the  lift-off  height. 

Until  this  point,  the  function  b2  was  used  to  extract  information 
on  the  ignition  and  growth  statistics.  The  reverse  direction  is  now 
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taken,  showing  how  h2  and  hi  depend  on  those  statistics.  Starting 
from  Eq.  (10),  and  taking  the  inverse  of  the  gradient  of  Pb2  at  h2 
proportional  to  L1/2,  the  following  expression  is  obtained: 


Zb, mm  +  k 2  [0  (fi gn;2> 


-1/2 


with  k2  a  non-dimensional  prefactor.  To  derive  an  expression  for  hu 
not  the  expected  number  of  ignition  events  per  time  per  axial 
length  is  needed,  but  those  per  time  per  surface  area  (presumably 
the  iso-surface  of  the  most  reactive  mixture  fraction),  because  the 
flame  propagates  on  a  surface  (the  iso-surface  of  stoichiometry). 
This  frequency  density  per  surface  area  /ign;l  scales  as  /ign;2  r~\  with 
r  the  characteristic  radius  of  the  axisymmetric  flame  surface.  As  the 
convolution  is  now  performed  over  a  surface  and  over  time  (instead 
of  over  a  length  and  over  time,  as  in  Eq.  (6)),  the  relation  between  hi 
and  the  ignition  and  transport  properties  becomes  more  complex. 
An  additional  non-dimensional  parameter  should  be  introduced, 
for  instance  On  dimensional  grounds  we  therefore  expect: 

‘'back 


hi 


Zb;min  +  ki[&  (fign;l>] 


with  /<i  a  non-dimensional  prefactor.  Note  that,  as  the  values  of 
Pm  (z)  can  never  exceed  those  of  Pb2(z),  hi  must  always  exceed  h2. 
This  imposes  limits  on  the  validity  of  Eq.  (14),  related  to  the  ratio 
of  the  axial  height  of  the  considered  region  and  the  circumference 
of  the  considered  flame  surface.  Solving  the  exact  equation  of  Pm 
for  pockets  growing  with  fixed  speeds  on  a  cylindrical  surface 
shows  that  the  function  is  strictly  decreasing  in  the  part 

\  ‘'back  / 

of  its  domain  exceeding  one.  The  trends  of  hi  and  h2  as  a  function 
of  jet  Reynolds  number  -  initially,  both  decrease  by  similar  degree, 
whereas  at  higher  jet  Reynolds  numbers  hi  increases  much  stronger 
than  h2-can  now  be  understood  through  these  arguments,  as  at  low 
jet  Reynolds  number  mainly  zb.min  is  affected,  whereas  at  higher 
Reynolds  numbers  the  different  scalings  cause  hi  and  h2  to  diverge. 


4.  Conclusions 

By  analysing  high-speed  recordings  of  luminescence  in  the  visi¬ 
ble  and  near  UV-part  of  the  spectrum,  the  lift-off  behaviour  of  flames 
burning  in  a  hot  and  diluted  coflow  was  studied.  Analysis  of  the  lumi¬ 
nescence  images  showed  that  the  physical  mechanisms  governing 
the  lift-off  process  in  jet-in-hot-coflow  flames  and  conventional 
lifted  flames  are  very  different.  This  difference  is  reflected  in  the 
trends  for  lift-off  height  as  a  function  of  jet  Reynolds  number,  which 
are  completely  dissimilar.  Instead  of  flame  propagation,  ignition  ker¬ 
nel  generation  by  autoignition  followed  by  convection  and  growth 
are  responsible  for  flame  stabilisation.  In  order  to  quantify  the  find¬ 
ings,  the  ignition  process  was  reduced  to  two  parameters:  the  axial 
location  where  ignitions  first  occur  zb;min  and  a  mean  ignition  fre¬ 
quency  density  in  the  stabilisation  region  </ign;2).  The  axial  flame 
growth  could  be  described  by  a  single  parameter  Based  on  obser¬ 
vations  and  dimensional  reasoning,  the  lift-off  height  follows  a  scal¬ 
ing  rule  as  described  by  Eq.  (14). 

Addition  of  higher  alkanes  and  increasing  the  coflow  temperature 
have  a  similar  effect.  Both  reduce  the  chemical  timescale,  leading  to  a 
lower  zb,min  and  higher  ignition  frequencies  (fjgn;2),  reducing  lift-off 
height.  The  influence  of  fuel  composition  or  coflow  temperature  on 
flame  pocket  speeds  is  relatively  small.  The  effects  of  jet  Reynolds 
number  on  the  parameters  are  more  intricate.  An  increase  in  Rey¬ 
nolds  number  from  3000  to  5000  lowers  in  all  cases  the  location  of 
first  ignition  zb,min.  This  is  likely  related  to  the  faster  entrainment  of 
the  hotter  part  of  the  coflow  for  increasing  Reynolds  numbers.  An  in¬ 
crease  in  jet  Reynolds  number  leads  to  significantly  higher  mean 
strains  of  the  axial  velocity  component  at  the  location  of  the  reaction 
zone,  as  shown  by  the  decrease  in  the  Taylor  length  scale.  At  jet  Rey¬ 
nolds  numbers  exceeding  5000,  this  is  accompanied  by  diminishing 


ignition  frequencies  in  the  flames  burning  in  the  colder  coflow,  while 
in  all  flames  the  growth  of  newly  formed  ignition  kernels  is  hindered. 
Extinction  of  ignition  kernels  is  occasionally  witnessed  at  the  highest 
Reynolds  numbers.  A  trend  seen  in  both  zb.min  and  (fign;2),  and  conse¬ 
quently  in  the  lift-off  height,  versus  the  jet  Reynolds  number  is  that 
increased  levels  of  turbulence  hinder  the  autoignition  process  more 
strongly  in  flames  with  slower  chemistry. 


Appendix  A.  Supplementary  data 

Supplementary  data  associated  with  this  article  can  be  found,  in 
the  online  version,  at  doi:  10.1 01 6/j.combustflame.201 0.01. 002. 
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